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CZ Abstract 

^^ This paper is devoted to the study of a fractional version of non- 

I hnear ^"{t), t > 0, hnear M^it), t > and sublinear 931" (i), t > 

death processes. Fractionality is introduced by replacing the usual integer- 
order derivative in the difference-differential equations governing the state 

I I probabilities, with the fractional derivative understood in the sense of 

Q^ Dzhrbashyan-Caputo. We derive explicitly the state probabilities of the 

Q, three death processes and examine the related probability generating func- 

• tions and mean values. A useful subordination relation is also proved, al- 

'"Ti lowing us to express the death processes as compositions of their classical 

C^ counterparts with the random time process T2v{t), t > 0. This random 

C time has one-dimensional distribution which is the folded solution to a 

' ' Cauchy problem of the fractional diffusion equation. 
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We assume that we have a population of uq individuals or objects. The compo- 
nents of this population might be the set of healthy people during an epidemic 
or the set of items being sold in a store, or even, say, melting ice pack blocks. 
However even a coalescence of particles can be treated in this same manner, 
^^ leading to a large ensemble of physical analogues suited to the method. The 

T-H main interest is to model the fading process of these objects and, in particular, 

h^ to analyse how the size of the population decreases. 

• rH The classical death process is a model describing this type of phenomena 

/\^ and, its linear version is analysed in fl], page 90. The most interesting feature 

of the extinguishing population is the probability distribution 



^ 



Pk{t)^ Pi- {M{t) = k\M{0)^no}, t>0,Q<k<no, (1.1) 

where M{t), t > is the point process representing the size of the population at 
time t. If the death rates are proportional to the population size, the process is 



called linear and the probabilities ( 1.1 ) are solutions to the initial- value problem 

' liPki*) = Kk + l)Pk+iit) - fikpk{t), < fc < no, 

,n^ jl' fc = "o, (1-2) 

0, < K < no. 



withp„o+i(i) = 0. 



The distribution satisfying (1.2) is 



Pkit) 



'fikt 
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-fit\^o-k 



< fc < no. 



(1.3) 



The equations ( 1.2 ) are based on the fact that the death rate of each component 



of the population is proportional to the number of existing individuals. 

In the non-linear case, where the death rates are /ifc, < fc < ng, equations 



(1.2 1 must be replaced by 



dtPk[t) = Mfc+iPfc+i(<) - fJ-kPhit), < k <no, 
1, fc = no, 
0, < fc < no. 



Pk{0) 



(1.4) 



In this paper we consider fractional versions of the processes described above, 
where fractionality is obtained by substitution of the integer-order derivatives 
appearing in (1.2 1 and (1.4 1, with the fractional derivative called Caputo or 



Dzhrbashyan-Caputo derivative, defined as follows 



d"f{t) _ 1 



/'(^) 



it 



ds. 



J'{t). 






(1.5) 



The main advantage of the Dzhrbashyan-Caputo fractional derivative over the 
usual Riemann-Liouville fractional derivatives is that the former requires only 
integer-order derivatives in the initial conditions. 

The fractional derivative operator is vastly present in the physical and math- 
ematical literature. It appears for example in generalisations of diffusion-type 
differential equations (see [2], [s], [2] and [5]), hyperbolic equations such as tele- 
graph equation (see [6|), reaction-diffusion equations (see [?]), or in the study 
of continuous time random walks (CTRW) scaling limits (see [s], lo]). Frac- 
tional calculus has also been considered by some authors to describe cahotic 



Hamiltonian dynamics in low dimensional systems (see e.g. [10], [TI], [12], 13 
and 14 ). For a complete review of fractional kinetics the reader can consult 
15 or the book by Zaslavsky 16 . In the literature are also present fractional 
generalisations of point processes, such as the Poisson process (see [l7] 
and 



, [20], 21 and [221) and the birth and birth-death processes (see [23], 
). Fractional models are also used in other fields, for example finance ( 

)■ 

The population size is governed by 
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^Pk(t) = ^J.k+lPk+l{t) - ^ikPkit), 



< fc < no, 



pm 



fc = "-0, 
< fc < no. 



(1.6) 



and is denoted by 9A'^{t), t>Q. 



Let us assume that a crack has the form of a process T2u{t), t > 0. For v = 
1/2, this coincides with a reflecting Brownian motion and has been described 



and derived in 28 . For v ^ 1/2, the process T2u{t), < > 0, can be identified 
with a stable process (see for details on this point [29]). The ensemble of no 
particles moves on the fracture and, at the same time, undergoes a decaying 
process which respects the same probabilistic rules of the usual death process. 
For the number of existing particles, we have therefore 

/•OO 

KW- / Pk{s)Pr{T2At)^ds}. (1.7) 

Jo 

We observe that 

Pi {T2At)eds}^q{.3,t)ds, (1.8) 

is a solution to 

Q2^ Q2 

qis,t) = -^q{s,t), s>0, i>0, (1.9) 



with the necessary initial conditions. Furthermore we recall that 

e-''^q{s,t)dt^z''-^e-''"', z > 0, s > 0. (1.10) 



The distribution q{s,t) is also a solution to 

■^qis,t) = -—q{s,t), s>0, (1.11) 



as can be ascertained directly. If we take the fractional derivative in (1.7) we 
get 



di 



J /"DO piV 

^Pk{t)=J Pkis)^PT{T2At)eds} (1.12) 

d 

Pk{s)^q{s,t)ds 

-q(s,t)Pk(s)\„ + / ^^ q(s,t)ds 



[-^J■kPk{s) + Hk+iPk+iis)] g(s, t)ds 

= -^ikPkit) + f^k+iPk+i{t)- 

This shows that replacing the time derivative with the fractional derivative 
corresponds to considering a death process (annihilating process) on particles 
displacing on a crack. 

We now give some details about ( |1.11[ ). By taking the Laplace transform of 
both members of (|1.11| we have that 



^td" ^^.,,,,_ d 



■OO 

e~^'^q{s,t)dt^-^{J e-^'q{s,t)dt] (1.13) 



(z^-'^e-''-'^) 



^2u-l^-sz-^ 



Furthermore, 



at-^ 



q{s,t)dt^z'' / e~^*q(s,t)di-z''"^q(s,0) 



(1.14) 



z" (z'^-ie-"^") - z'^'His), 



and therefore for s > this estabhshes that q{s,t) solves equation (1.11). We 



note that a gas particle moving on a fracture has inspired to different authors 
the iterated Brownian motion (see pJO]). 
The distribution 



pl{t) = Pr {M^it) = fc I 5M"^(0) = no} , < A: < ng, 

is obtained explicitly and reads 



n ^^3 2^ -^ ' < fc < no, 



(1.15) 



pUt) = { 



h=k 
no Tio 

^ E^A^f-^rnt"), fc = 0, no > 1- 



(1.16) 
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m—l h—1 

h^rn. 



M/i - ^^r. 



Obviously, for /c = 0, ng = 1, 



The Mittag-LefFler functions appearing in (1.16) are defined as 



h=0 



r {vh + 7) ' 



J^,7 > 0. 



(1.17) 



(1.18) 



For V = ^ = 1, _Bi^i (cc) = e^ and formulae ( 1.16 ) provide the explicit distribution 
of the classical non-linear death process. 

For /ife = kjjL the distribution of the fractional linear death process can be 



obtained either directly by solving the Cauchy problem (1.6) with /i^ — k ■ fi 



and Pno+i(t) = 0, or by specialising (1.16) resulting in the following form 

pUt)- (^l^yf: ("°;''){-irE^A-ik + r)K). (1.19) 

A technical tool necessary for our manipulations is the Laplace transform of 
Mittag-Lefflcr functions which we write here for the sake of completeness: 



e-'H-'-^E,,^{±m-')dt = -, «n(z) > 1^1^. (1.20) 

z'^ ^v 



Another special case is the so-called fractional sublinear death process (for 
sublinear birth processes consult Isl]) where the death rates have the form 
jUfe = /x(no + 1 — fc). In the sublinear process, the annihilation of particles or 
individuals accelerates with decreasing population size. 

The distribution p^(i), < A; < no of the fractional sublinear death process 
S[Tt'^(i), t > 0, is strictly related to that of the fractional linear birth process 
N^it), t > (see, for details on this point, [24]): 

Pr{9Jt''(i)=0|97l''(0)=no} = Pr{iV''(i) > no | N''{0) = 1}. (1.21) 

In general, the connection between the fractional sublinear death process and 
the fractional linear birth process is expressed by the relation 

Pr {M'^it) == no - (fc - 1) I W(0) = no} (1.22) 

= Pr{iV'^(i) = fc I Ar'^(O) = 1}, l<fc<no. 

This shows a sort of symmetry in the evolution of fractional linear birth and 
fractional sublinear death processes. 

For all fractional processes considered in this paper, a subordination rela- 
tionship holds. In particular, for the fractional linear death process we can write 
that 

Ar{t) = M{T2y{t)), < !^ < 1, i > 0, (1.23) 

where T2^{t) is a process for which 

Pr {T2^(t) e ds} = q{s, t)ds, (1.24) 

is a solution to the following Cauchy problem (see [32] ) 

^g(s, t) = ^9(s, t), t > 0, s > 0, 

i9(s,i)L=o = 0' (1-25) 

j{s,0)^S{s), 0<zy<l, 

with the additional initial condition 

gt(s,0)=0, 1/2<J.<1. (1.26) 



In equation (1.23), M{t), t > 0, represents the classical linear death process. 



Subordination relations of this type are extensively treated in 29 1 and 33 . 

We also show that all the fractional death processes considered below can 
be viewed as classical death processes with rate fi ■ S, where S is a Wright- 
distributed random variable. 



2 The fractional linear death process and its 
properties 

In this section we derive the distribution of the fractional linear death process 
as well as some interesting related properties and interpretations. 
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Figure 1: Plot of p"-^(i) (in black) and p^ (t) (in grey), both with no = 10. 



Theorem 1. The distribution of the fractional linear death process M'^{t), i > 
with riQ initial individuals and death rates /x^ = /i • fc, is given by 



Pl{t) = Pr{Ar{t) = k I Ar{0) = no} 

'''^^^f:h;%-irE.A-ik+r),tn. 



(2.1) 



where < fc < no, i > and v G (0,1]. The function E^^i{x) is the Mittag- 



Lejfler function previously defined in (1.18) 



Proof. The state probability PnAt), t > is readily obtained by applying the 
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Laplace transform to equation ( |1.6[ ), with /i/j = H ■ k, and then transforming 
back the result, thus yielding 



K„(t) = E^^^i-n^tit-), t > 0, z. e (0,1]. 



(2.2) 



When fc = no — 1, in order to solve the related differential equation, we can 
write 



v!^-l 



Z^CiPno-l) [z) = M^O-^T— M("-0 - l)^{Pno-l} (^) 



z" + no^ 



(2.3) 



<^ C {pno-i} {z) ^ tiUaz" 



<(=> £{p„o_i}(z) = noz'' 



1 



1 



z" + no^ z^ + (no — 1)^ 
1 1 



z" + (no + l)/i z" + no/i 




Figure 2: Plot of p°-J_j(i) (in black) and p^g^iit) (in grey). Here no = 10. 



By inverting equation ( |2.3[ ), we readily obtain that 

Ko-i W = no (-B.,i(-(no - 1)^^^") - S.,i(-'^oAit")) • (2.4) 



For general values of k, with < /c < no, we must solve the following Cauchy 
problem: 



|,p.w^M^^+i)(;;^ 



(2^5) 



no — k — 1 , I, 1 ^ 

E l °" " )(-l)'^i?.,i(-(fc+l + r)Ait"^)-/ifcpfc(0, 



r=0 



subject to the initial condition pa:(0) — and with v e (0, 1]. The solution can 
be found by resorting to the Laplace transform, as we see in the following. 



z''/:{Pfc}(z) = A*(fc + i) 



no 
/c + 1 



rio-fc-l / 7 1 > 



r=0 



I/-1 



(2.6) 
^fc£{pfe}(z). 



z"" + (fc + 1 + r)^ 
The Laplace transform £{pfe} (z) can thus be written as 
C{Vk}{z) (2.7) 



"•'-"UTi S 



no \ ■r-^ /no — K — 1 



r=0 



(-1) 



vl'-l 



1 



riQ — fc — 1 / , 

no \ v-^ ( uq — k 



E 



k J ^—' \r^\} Kz'' + kn z" + (k + 1 + r)u, 



z" + {k + l + r)fi z"" + k^ 
1 1 



(:)f("°:>-'r'-'(;4 



1 



7iQ — k 



+ k^i z'' + {k+j)^i 



^u-l 






no\ z ■^-^ /no — fc 



k J z'^ + kii 



E 



J 



(-IF 



UQ — k 



^u-l 



nQ\ i-^ I no - fc\ z" 

^^ V i / z- + (fc+j)/i- 



no\ z 



u-l 



no\ y- /»0-fcV ^y, ^ 

^ \ j / z'' + (k + j)fi ' \ k J z" + kfi 



n 



5?/ taking now the inverse Laplace transform of (2.7), we obtain the claimed 
result (ll.igp. 



Remark 1. When v = 1, equation (1.191 easily reduces to the distribution of 
the classical linear death process, i.e. 



Pk{t) 



'^Me-''^^*(l-e-'^*)"""' 



t > 0, < fc < no. (2.i 



In the following theorem we give a proof of an interesting subordination 
relation. 

Theorem 2. The fractional linear death process M'^(t), t > can be represented 
as 

Ar{t) '^- M{T2,{t)), i > 0, z. e (0, 1], (2.9) 

where M{t), t > is the classical linear death process (see e.g. W, page 90) 
and T2i/{t), t > 0, is a random time process whose one- dimensional distribution 
coincides with the folded solution to the following fractional diffusion equation 



' ^q{s,t)^£,q{s,t), t>0, :.e(0,l], 
q{s,0)^5{s), 



(2.10) 



with the additional condition qt{s,0) — if ly E (1/2, 1] (see [32|j. 



Proof. By evaluating the Laplace transform of the generating function of the 
fractional linear death process M^{t), i > 0, we obtain that 



e-^'G''{u,t)dt 



(2.11) 






^^-^ \ k / ^^ \ r 



r!^-l 



/"CXD ^0 / \ UQ — k / , \ 



fe=0 



>fe=o 

"0 






^^S^^ ^V—\ 



z^ G[u^ s)ds 



/ ''^'l ^u''Pr{M(s)^k}fT,As,t)dsdt 
Jo Jo ^^f) 

poo ( °° "I 

J '^^'' \ E "'P^ {^^(^2.(t)) = fc} I di, 



and this is sufficient to prove that (2.9) holds. Note that we used two facts. The 
first one is that 



e~"*/T2„(s,0^^ = 2''"^e- 



s > 0, z > 0, 



(2.12) 



is the Laplace transform of the solution to (2.101. The second fact is that the 



Laplace transform of the Mittag-Leffier function is 



1° 

Jo 



e-^*£;^,i(-^r)di 



z'^ + t?' 



(2.13) 



n 



In figures [T] and [2J we compare the behaviour of the fractional probabilities 
PnoW ^^'^ PnJ-ii't) ■^ith their classical counterparts Pn^{t) and Pn^-iit), t > 0. 



What emerges from the inspection of both figures is that, for large values of t, 
the probabilities, in the fractional case, decrease more slowly than p^^ (i) and 
Png-i{t). The probability p°J_i(i), increases initially faster than p^^_;^(i), but 
after a certain time lapse, p,\ _i(0 dominates p^j^„]^(t). 

Remark 2. For v — 1/2, in view of the integral representation 



EiAx) 



— u; -^2xw 



dw. 



X G 



we extract from (|1.19| that 

-pUy)dy^Pr{M{\Bm^k}, 



g-^^(^o)^ Y^ (^"0-fc'||-_-^^rg-2to(fe+r)MtS 



E 

r=0 



r 



(2.14) 



(2.15) 



where Bit), t > is a Brownian motion with volatility equal to 2. 



Remark 3. We can interpret formula (1.191 in an alternative way, as follows. 

(2.16) 



For each integer k G [0, no] '"'fi have that 

pl{t) = Pr{IvF{t) = k\Nr{0) = n^} 

= / Pk{s)Pr{T2.{t) e ds} 
Jo 

= (T) "E ("°; y^iyfe'^C^^^^PrlnAt) G ds} 

= (T) 'e (^°;''){-ir l^ e-^^'^'-^H-^w^^^^M-st-nds 

W^-.a-.(-0^^{A^s(i") = k \ M^{0) = no}de, 



where W-v.i-u{^C) ^^ ^^^ Wright function defined as 

{-if 



w^-..i-.(-o-E 



^_.Hr(l-t.(r + l))^ 



Q <v<l. 



(2.17) 



FFe therefore obtain an interpretation in terms of a classical linear death process 
Msit), i > evaluated on a new time scale and with random rate fj, ■ S, where 
'R is a random variable, ^ G M"*", with Wright density 



/s(0 = M/-.4-.("0, 



ee 



(2.18) 
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From equation (1.6) with fik ~ k ■ /i, the related fractional differential equa- 



tion governing the probability generating function, can be easily obtained, lead- 
ing to 



^" G-'iu, t) - -^Jiu{u - l)§-^G-{u,t), V e (0, 1], 



dt 



G"'(u,0) = u" 



(2.19) 



From this, and by considering that E^f''(i) = ■^G'^iu^t')^ , we obtain that 



^EM''(t) = -[iKWit], V e (0, 1], 
EM''(i) == TiQ. 



(2.20) 



Equation (2.20) is easily solved by means of the Laplace transforms, yielding 
EM''(t) = noE^.^ii-nf), i > 0, i/ e (0, 1]. 



Remark 4. The mean value EM'^{t) can also be directly calculated. 

"0 

EArit) = J2kpUt) 

> "0 / 7 \ 



(2.21) 



(2.22) 



fc=0 
no 



fc=0 
no 



EM. E 



r ~ k 



i-ir-^E^^ii-rfitn 



j2E.A-rK)i-irJ2'' 



r=0 

no 



fc=l 



no \ / no 



r — k 



{-If 



r=l ^ '' k=l ^ ^ 

= noE^^i{-^it''). 



This last step in (2.22) holds because 

r-l 

r 
k 



fe=l ^ ^ /c=0 



10, r>l. 



(2.23) 



3 Related models 

In this section we present two models which are related to the fractional linear 
death process. The first one is its natural generalisation to the non-linear case 
i.e. we consider death rates in the form fik > 0, < k < uq. The second one 



is a sublinear process (see 31 ), namely with death rates in the form /i^ 



yu(no -|- 1 — A;); the death rates are thus an increasing sequence as the number of 
individuals in the population decreases towards zero. 



11 



3.1 Generalisation to the non-linear case 

Let us denote by fW*^ (i) , t > the random number of components of a non-linear 
fractional death process with death rates fJ-k > 0, < k < uq. 

The state probabilities p^(i) = Pr {fM^it) ^ k\ fM^^O) ^ no}, i > 0, < fc < 
no, v G (0, 1] are governed by the following difference-differential equations 

^pfe(t) = ^j.k+iPk+i{t) - ^J-kPkit), < fc < no, 



dt 

di'^PnoV') = ^t^noPnoV'jj 



k = 0, 

k = no. 



(3.1) 



PkiO) = 



0, < fc < no, 

1, k — uq. 



The fractional derivatives appearing in (3.11 provide the system with a global 



memory; i.e. the evolution of the state probabilities Pkit), i > 0, is influenced 
by the past, as definition (1.51 shows. This is a major difference with the 



classical non-linear (and, of course, linear and sublinear) death processes, and 
reverberates in the slowly decaying structure of probabilities extracted from 



(3.1) 



In the non-linear process, the dependence of death rates from the size of 
the population is arbitrary, and this explains the complicated structure of the 
probabilities obtained. Further generalisation can be considered by assuming 
that the death rates depend on t (non- homogeneous, non-linear death process). 

We outline here the evaluation of the probabilities p'^{t), t > 0, < k < Uq, 
which can be obtained, as in the linear case, by means of a recursive procedure 
(similar to that implemented in f24' for the fractional linear birth process). 



Let k — Uq. By means of the Laplace transform applied to equation (3.11 
we immediately obtain that 



flW=^..l(-/^noi'')- 



(3.2) 



When fc = no — 1 we get 






<^ ^ {Pno-l} i^) ^ f^no 



1 






Z'' + flng Z" + ^no-1 

1 1 



Z -\- fJ^iiQ ^ ~r ^riQ — l 



1 



(3.3) 



A^no — 1 Mno 



Eu,l{~fJ'not'') ~ E^^{~^no-lt^) 



For fc = no — 2 we obtain in the same way that 



(3.4) 
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= -t^no-2^{Pno-2}iz) + 



MnoA^no — 1 



SO that 






(3.5) 



A^no — 1 A*no 

Mno — 1 ^ A*no 
1 



1 



1 



Z -f- /iriQ Z + f^riQ — l 

1 1 



1 



1 



Z -\- I-Luq ^ ~r /^no — 2 / Mno— 2 Mno 

1 \ 1 



^z'^+/i„Q_l Z^ + lJ"na-2 J A'no-2 ^ A^rio-l 

By inverting the Laplace transform we readily arrive at the following result 



Pno~2W — MnoMno 



(Mmo-I " ^^no)i^^no~2 - Mno ) 



(3.6) 



(/-^no-l ^ ^J■no){^^no-2 - A^rio ) (A^mo-I ^ A^no)(A'no-2 ^ A^no-l) 

-E',y,l( — /i„(,-2^'') 



-£'i>,l(^Mno^'') 



A'noMno — 1 



(A^no-1 — /^rio)(Mno-2 ^ Mno ) 



(.Mrao — 1 Mno/ VMrio— 2 MriQ — 1 Mno — 2 A^no 



(/^rio-1 — Mno)(Mno-2 ^ MriQ-l) 



A'noA^no — 1 



(.A'no — 1 A*no/vA^no — 2 Mno / 
-E'i',l( — Mno-2^'') E,y^i{ — IIn„-lt'') 



(A«no-2 - Mno-l)(A^no-2 - A'riQ ) (Mmo-I ^ ^^no)i^^no-2 - A«no-l) 

The structure of the state probabilities for arbitrary values of fc = ng — /, < 
I < no, can now be easily obtained. The proof follows the lines of the derivation 
of the state probabilities for the fractional non-linear pure birth process adopted 
in Theorem 2.1 in f24l. We have that 



(-1 



E^^li^lJ.no-mt'^) 



pL-i{t) = <^ 



11 Mno-j ^/^ ; 

~ 11 (f^no — h^^^no — m) 



h=0 



.Ei.^li — finot"), 



I <l <no, 



1^0. 



(3.7) 



13 



By means of some changes of indices, formula (3.7) can also be written as 



Pkit) = <^ 



11 f^J 2^ no ' 

^='= + 1 ra=k n {^Jih-^^n^) 
h=k 


< k <no, 


[E,A-^^non, 


k = Uq. 



(3.8) 



For the extinction probability, we have to solve the following initial value 
problem: 



^po(t) = Ml ri Mi Yl 



^ ^ E^^.i-flrntn 



na 



h=l 

^po{t) = fiiE,^i{~fiin, 

Ifo(O) - 0, 



no > 1, 



no = 1, 
no > 1. 



(3.9) 



When no > 1, starting from (3.91 and by resorting to the Laplace transform 
once again, we have that 



"0 "" 1 r-l 

c{fo} {z) = n M, E -. ^7X77 

h=l 



(3.10) 



The inverse Laplace transform of (3.10) leads to 



poit) - n H E —0 

■J 1 m=l J] (^f^ - n^) 



h=l 



1 



i'^ E^^^j^l(~ flrnf^) 



1 



(3.11) 



= n M, E -1^0 — 7^ [1 - E.^{-^^ran 

^ = 1 m=l n (M/.-Mm) ^"^ 
h=l 

h^rn 



En 



M/7 



m— 1 /i^l ^' ' ^ rn— 1 h^l ^' 



no no 



En 

hy^ni 



^J■h 



£'i/,l( — Mmi") 



1-En 

m— 1 h—1 



fJ-h 



IJ-h - fJ-r. 



Ei^^li — flmf^)- 



Note that, in the last step, we used the following fact: 

no no / 

771—1 h—1 ^ 

hy^ni 



(3.12) 
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This can be ascertained by observing that 



where 



[ ifj-h 


- w) = 


dotA = 


= > >i,.(-i)' 


+iMini 


l<h<l<no 






i=i 




1 


1 


1 




A = 




M2 


• • • M^o 


, 




tH"-' 


i^T-' 


„no-l 
Mno 





(3.13) 



(3.14) 



is a Vandermonde matrix and Mini j is the determinant of the matrix resulting 
from A by removing the first row and the j-th column. 
When no = 1 we obtain 



so that the inverse Laplace transform can be written as 

;^^(t)=Aiii"^^,^+i(-Mii") 
= 1-S,,i(-Mit"). 

We can therefore summarise the results obtained as follows: 






■no ,^_:^ 

n M,E 



< k < no, no > I, 



k = no, no > 1, 



and 



fJ-h 



"0 "0 

1 - E n ( -^^^^- ) E^.ii-l^^mn, no > 1, 



.l-^..i(-A*in, 



no = 1. 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



3.2 A fractional sublinear death process 

We consider in this section the process where the infinitesimal death probabili- 
ties have the form 



Pr {m{t, t + dt] = -1\ mit) =k} = fi{no + 1 - k)dt + o{dt), 



(3.19) 



where ng is the initial number of individuals in the population. The state 
probabilities 



Pfc(t) = Pr{9Jl(i) = k I »t(0) = no} , < fc < no, 



(3.20) 
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Time (arb. units) 



Figure 3: Plot oi p^J{t) (in black) and p„J(i) (in grey), with no = 2 



satisfy the equations 

^pfe(i) = -Ai(no + 1 - k)pkit) + ^(no - fc)pfe+i(t), 1 < fc < Hq, 



''po(0 = M»^opi(t), 



dt 

Pfc(O) 



0, 



|l, k^no, 
lO, 0<fc<no. 



(3.21) 



Pfc(O) = 



(3.22) 



In this model the death rate increases with decreasing population size. 

The probabilities p^(t) = Pr {m^it) = k \ 97l''(0) == no} of the fractional ver- 
sion of this process are governed by the equations 

' w^Pk{t) = -Ai("-o + 1 - k)pk{t) + A*("o - k)pk+i{t), I <k <no, 
§7po{t) ~ M'^opi(i), fc = 0, 

1, k ~ no, 
^0, 0<fc<no. 

We first observe that the solution to the Cauchy problem 
fl^P«o(i) = -MPnoW> 

In order to solve the equation 
|p„„_i(0)-0, 



(3.23) 



(3.24) 
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we resort to the Laplace transform and obtain that 



z" + n z"^ + 2/i 
1 1 



z"" + ^1 z'^ + 2^ 



By invertmg (3.251 we extract the foUowmg result 

Pno-iW = E,,i{-K) - i?.,i(-2Mt''). 
By the same technique we solve 

f £pn„-2{t) = -3fiPno-2{t) + 2/i [E,,i{~ fit") - ^,,i(-2/ir)] ^ 

lp„„_2(0)-0, 



thus obtaining 



^{p:;„-2}w = 2/.z''-i 
1 



1 



1 






z'' + /i z"^ + 2/i 
1 \ 1 



1 



z^ + zi z'' + 3/i/2/i Vz'^ + 2/i z^ + S^Jfj. 



z" + 3/i 
1 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



1 \ 1 



,i/-i 



,y-i 



-2- 



2"^+/* 2"^+ 2/i Z-^+S/i 

In light of ( 3.28[ ), we infer that 

p«o-2W = KA-^^n - 2s.,i(-2/it'') + E,^,{-3^Ie) 



(3.29) 



For all 1 < riQ — ?7i < no, by similar calculations, we arrive at the general 
result 

P,VmW=E 7 (-l)'^-.i(-(^ + l)'^*'^)' l<"o-m<no. (3.30) 



i=o 



Introducing the notation uq ^ m ~ k, we rewrite the state probabilities (3.301 
in the following manner 



no — fc 



put) =T.(''\ '') (-l)'^-i (- a + 1) f^n , l<k<no. (3.31) 



!=0 



/ 



For the extinction probability we must solve the following Cauchy problem 

(3.32) 



'^■' po(t) = M"o Er=o'' (T') (-1)'^-,! (- (l + 1) Mt"-) , 



dt 

Po(0)=0 
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Time (arb. units) 



Figure 4: Plot of p"J_^(t) (in black) and p"'7_i(t) (in grey), with no = 2. 



The Laplace transform of (3.321 yields 

"0-1 



z''C{p'^,}(z)^^in,Y. 



1=0 



^u-l 






(3.33) 



The inverse Laplace transform can be written down as 

"0 — 1 / -1 \ 1 nt 



1=0 



The integral appearing in (3.34) can be suitably evaluated as follows 



(3.34) 
(3.35) 







E 

oo 

E 



(-a + iv) 



m rt 



s^'^it-sf-'ds 



T{vm + 1) Jq 

{-{I + 1)^*)" i'-'("+i)r(i/)r(i'm + 1) 

T{i>'m + 1) T{i>m + 1^+1) 



r(^) 



D) ^. 



{-{i + i)fit")"^+^ 



i-Ki + ^))^, r(Km + i) + i) 
^-J>l^K.(-(/ + i),n-i] 
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By inserting result (3.35) into (3.34), we obtain 



"0-1 



I 

no-] 

= E 



1=0 
no-1 



I 



l + l 



1=0 

no 



l + l 






(3.36) 



1=1 



E 7 (-i)'^^^^H/^^^)-E ; M) 



1=1 



1=1 



i+j:r){-iyE.A-i^^n 



1=0 



E(7)(-i)'^-iH^^^)- 



Remark 5. We check that the probabilities (3.311 and (3.361 sum up to unity. 
We start by analysing the following sum: 



no no no — fc ^ ,\ 



k=l 



k=l 1=0 



I 



In order to evaluate (3.37), we resort to the Laplace transform 

jj_l no no-fe 

— ^ "C 

I 



E/:{P^}w = VEEf"°r'V-i)'T ' 



k=l 



M 



k = l 1=0 



l + l 



(3.37) 



(3.38) 



By using formula (6) of [34| (see also 35 , formula (5.4-1), page 188), we obtain 
that 



"0 y_l "0 

z — 



r(no-fc + l) 



^ ^ tt(^ + l)(f +2)...(^ + l + no-fc) 



(3.39) 



k=l 



,.-i^T[^ + l)V{n,-k + l) 



M 



fc=i r(^ + l + no-A:) 
E / a;^(l-a;)"»"'=da; 



M 



1,_1 "0 ^1 



fe=l 



= / xV-i[l-(l-x)"«]dx 

z t^ Jo 

z M Jo 



a;7r-i(i_a;)"Vx 



e ^ f' (l — e ^) " dy 
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{y/ll=w) 1 ^_i 

Z 



^ l,_n V ^ / ^0 



i:=0 
no 



A^ V i-y^ ^ z" + fik 



k=0 



The inverse Laplace transform of (3.391 is therefore 

"0 



fc=l 



fc=0 
no 



Y.Pkit) = l-J2( ° (-l)'i?.,i(-MA:n 



(3.40) 



A,-l 



5: 7 (-if£;,i(-Mn. 



By putting (3.361 and (3.40) together, we conclude that 

no 



Y.put) = h 



(3.41) 



k=Q 



as it should be. 



Remark 6. We observe that, in the linear and sublinear death processes, the 
extinction probabilities coincide. This implies that although the state probabil- 
ities p^{t) and p\.{t) differ (see figures 3 and 4) for all 1 < k < hq, we have 
that 



Y,pi{t) = Y,Pi[t). 

k=l k=l 

This can be checked by performing the following sum 



Y.C{pl{t)}{z) 



(3.42) 



(3.43) 



k=l 



"0 / \ no-fe 


V r J^ z''+n{k + r) 






,-1 ^ /„o\ 


(no-fc)! 


(f+^)(f + fc + l)...(f+no) 




r(no-fc + i)r(^ + fc) 


r(f + no + l) 
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U-1 fl "0 , , 



Ai Jo 



fc=l 



(i„^)V-i^ ""U"o-'=(i_xfd,c 



1 Z'--! '•I 



^ Ai Jo 



x"°{l-x)^^^dx. 



This coincides with the fourth-to-last step of (3.391 and therefore we can con- 
clude that 



no 



fe = l 

3.2.1 Mean value 



fe=i 



E^J^w = -E ni-i)''E.A~i^ktn = Ep^w- 



(3.44) 



fe=i 



Theorem 3. Consider the fractional sublinear death process dJl'^it), t > de- 
fined above. The probability generating function &'^{u, t) — X^felo ^'^Pfc(^); t > 0, 
\u\ < 1, satisfies the following partial differential equation: 

|^©^(«, t) = ^iino + 1) (^ - l) [®'^(«, t) - p-^M + m(u - 1)|^®'^(", i)- 

(3.45) 
subject to the initial condition (5''(w,0) = u"" , for \u\ < 1, t > 0. 



Proof. Starting from (3.22), we obtain that 

"0 



d^ 

dF 



E"'p^(*) 



(3.46) 



fc=0 

no "0 — 1 

-/iE"'(«o + 1 - k)pl{t) + /i E "'K - fc)p^+i(i), 

fc=l fc=0 



so that 

§^'&''iu, i) - - Kno + 1) [e'-iu, t) - p;i(i)] + l^u^&^iu, t) (3.47) 

+ ^^"""^^^ [^''iu,t) - p[;(i)] - M^®^(«,i) 
u ou 

= A^(no + 1) Q - 1) [^''{u,t) - p;;(t)] + ^,{u - i)^e^{u,t). 

u 

Theorem 4. The mean number of individuals E9Jt'^(t), t > in the fractional 
sublinear death process, reads 

E9Jl^(i) = E 1 ? 1 ) (-l)'+'^..i(-Mi''). t > 0, z. e (0, 1]. (3.48) 



fc=i 



fc + 1 
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Proof. From ^A5\ and by considering that EdJl^it) = £«"(", OL^i> we di- 
rectly arrive at the foUowing initial value problem: 



\Em''{0)=no, 
which can be solved by resorting to the Laplace transform, as follows: 



(3.49) 



C{Em''{t)}{z) =no- 



z" - ^ 



k=l 



no 



'■0 

E 



1 



z'^ + Ijk z^ ^ ^ 

(3.50) 



z- - u -^^ V A: + 1 ' ^ ' 

fc— 1 



z z 

z'^ — fi z^ + Ilk 



In ( 3.50 ) , formula (3.36 ) must be considered. By inverting the Laplace transform 

(3.51) 



we obtain that 
EM^it) 



noE^^^ifitn + E f? ^!) (-1)' [E.Af^n - E^^.i-fiktn] 



fc=l 



fc + 1 



= noE,,,{ptn + E,s{K) E f? ^^') (-1)' 



k=l 



fc + 1 



±h^l\-^fE.A^^^kn 



k=l 



fc + 1 



^Ef?^'U-i)'^+^£^.i(-M^''), 



fc=l 



fc + 1 



as desired. 



n 



Remark 7. T/ie mean value (|3.48p can afeo 6e directly derived as follows. 

(3.52) 



fc=0 

"0 "Q- fc / 7 \ 

-E^E ("": )(-i)'i^M(-a+i)Mn 



no rio + l — fc 

E^ E 

fc=i ;=i 

no 



no - k 
l-l 



i-iy-'E^A-fiir) 



n p+l -l , , 

-E(-i)'"^-(-A'^^'') E h7-i 

Z=l fe = l ^ 
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It is now sufficient to show that 

na + l — l 



E ^ 



fe=i 



no — k 
l-l 



no + 1 



Indeed, 



no + l-' / 


-')- 


"0-1 

k=l-l 
"0 — 1 



EK-^)L , 



(3.53) 



(3.54) 



^ (no + 1 - fc - 1) 



fc=/-i 



A;=(-l 
no 



K + i)E UIi -'E 



"0-1 



k=l 



fc-1 



K+i)E ;:; -^E 



fe=i-i 

"0 + 1 



k=l + l 



k + l 

k-l 
I 



(».. 1)1 7 -<";;; 



no + 1 
/ + 1 



The crucial step of (3.54) is justified by the following formula 

riQ 



SG-0- 



3 
J-1 



no - 1 



(3.55) 



Figure [5] shows that in the subhnear case, the mean number of individuals 
in the population, decays more slowly than in the linear case, as expected. 

Note that ( |3.48 ) satisfies the initial condition EOT''(0) = no. In order to 
check this, it is sufficient to show that 



"0 + 1 



fc=l ^ ^ r=2 



no + 1 
r 



■"0 + 1 

E 

. r=0 



no + 1 



(-1)^ 



- 1 



(-ir 

no + 1^ 



(3.56) 



1 



= no- 



The details in (3.56) explain also the last step of (3.51) 



3.2.2 Comparison of dJV^{t) with the fractional linear death process 
M'^{t) and the fractional linear birth process N'^{t) 

The distributions of the fractional linear and sublinear processes examined above 
display a behaviour which is illustrated in Tablefll(see also Table 2 for the mean 
values) . 
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Time (arb . units) 



Figure 5: Plot of EM"'^(f) (in black) and E»t"'^(t) (in grey), no = 2. 



The most striking fact about the models dealt with above, is that the linear 
probabilities decay faster than the corresponding sublinear ones, for small values 
of k; whereas, for large values of k, the sublinear probabilities take over and the 
extinction probabilities in both cases coincide. The reader should also compare 
the state probabilities of the death models examined here with those of the 
fractional linear pure birth process (with birth rate A and one progenitor) . These 
read 



Pk 



3 = 1 ^■' ^ 



k>\. 



(3.57) 



Note that ji\{t) = E^^x{-W) is of the same form as '^"^^{t) = E^^ii-fxt"). We 
now show that 



k—nQ-\-l k—1 

no k 



fc=lj=l ^-^ ^ 

"0 "0 /J, _ 1 

i-E(-i)-^i?,,(-A,nEr:j 
j=i k=] ^-^ 

j=i V -^ / 



(3.36) with A replacing ii. 



(3.58) 
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Table 1: State probabilities p'^{t) for the fractional linear death process M'^{t), 
t > 0, and pfc(t) for the fractional sublinear death process DJl^{t). 

State Probabilities 

Ko(«) = ^..l(-M«'') 

Ko-i(*) = "0 [S^,i(-(no - 1)^") - E^,i{-nofit'')] 

Ko-iW = E.A-^^n - i?.,i(-2Mr) 

pm = (?) Er^-' ("v''")(~i)'^-,i(-{fc + OMt'o 
p%(t) = Er=v' rr')(-i)'-E.,i {- n + 1) /.t-) 

p5'(*) = "oEr=v' rr')(-i)'^-i(-(i+OMt'') 
pr{t) = Er^o"' ("v')(-i)'^-.i (- (' + 1) M") 

Po(t) = Er="o(7)(-i)'^-.iHMt'') 
Po{t) = Er=%(7)(-i)'^-iH/^*") 



Table 2: Mean values for the fractional linear birth N'^{t), fractional linear 
death M"{t) and fractional sublinear death dW{t) processes. 



Note that in the above step we used formula (3.55). 

By comparing formulae (3.4) of [24 and (3.31) above, we arrive at the con- 
clusion that (for \ = fi) 



PT{N''{t)^k\N''{0) = l} (3.59) 

-E(-:;)(-ir^^M(-A,n 

= Pr {Wl^it) = no + 1 - fc I 9^(0) = no} , 1 < fc < no. 



The probability of extinction po (t) corresponds to the probability of the event 
{N'^{t) > no} for the fractional linear birth process. 
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